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PRÓLOGO 


Alo largo de la historia los seres humanos siempre hemos aprovechado el agua como 
fuente para usos diversos y de todo tipo. La ingeniería hidráulica específicamente, 
ha servido para abastecer de agua y arbitrar las medidas de protección contra 
inundaciones, tormentas, procesos contaminantes, entre otros. 


En este tratado de ingeniería se puede encontrar los principios básicos de la hidráulica 
de canales abiertos y una guía académica de fácil comprensión para emprender el 
diseño de obras civiles de transporte de agua en lámina libre, como sistemas fluviales 
y agrícolas, redes de alcantarillados y otras aplicaciones. 


En la primera parte se encontrará un estudio de las ecuaciones fundamentales que 
gobiernan el flujo uniforme del agua. Las expresiones algebraicas se presentan con 
un breve desarrollo matemático, que fortalece el conocimiento del lector. 


De manera consecutiva, se presenta los métodos de cálculo para determinar 
el calado normal, así como los criterios de diseño de canales de contorno rígido 
y contorno erosionable, con base a secciones hidráulicas óptimas y velocidades 
máximas permisibles, respectivamente. 


Por último, se aborda el flujo no uniforme en canales hidráulicos mediante el 
desarrollo de conceptos como la energía específica de una sección, el tirante y 
pendiente crítica y el análisis de las superficies libres. 


El abordaje académico del trabajo se lo puede resumir en los siguientes tratados: 
flujo uniforme en canales, flujo no uniforme en canales y problemas de aplicación. Un 
resumen amplio y práctico para consolidar y mejorar el trabajo técnico profesional 
en esta rama de la ingeniería. 


La idea del autor es compartir sus experiencias y contribuir en la construcción de 
diseños de obras hidráulicas para el transporte de agua en lámina libre. 


Espero con modestia haber cumplido con mi cometido y afán académico. 


El autor 
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CAPÍTULO 1. FLUJO UNIFORME EN CANALES 


1.1. Introducción 


El fluido en canales abiertos es caracterizado por la presencia de la interface entre la 
superficie líquida y la atmósfera, tal es el caso del flujo en ríos, canales y alcantarillas, 
donde el flujo ocupa solo una parte de la sección. Por lo tanto, a diferencia del flujo 
en tuberías, la presión sobre la superficie del líquido en canales abiertos siempre es 
la presión atmosférica. 


Figura 1. Movimiento uniforme en cauces abiertos. 
Fuente: Elaboración propia. 


En la Figura 1 se presenta un fluido uniforme que corre por un canal abierto una 
longitud (1). La línea de gradiente hidráulico (L.G.H) se designa con la letra 1 y es 
igual a la relación entre la variación de energía entre las dos secciones y la distancia 
horizontal entre las mismas CL ), medida en el plano horizontal de referencia; 
mientras que, Í es la pendiente de la superficie libre y 3 es la pendiente de la solera 
del cauce. En canales con movimiento uniforme la línea L.G.H (gradiente hidráulico), 
pendiente de la superficie libre del cauce y la pendiente de la solera del canal son 
paralelas, es decir i= 1 = 1 . Esto se debe a que los parámetros: área (A ), velocidad 
(v) y gasto (Q ) son constantes en todas las secciones de la tubería. Es decir: 
i=I=1 


Ss 
dh, d(z+h) dz 
dl dl dl 
Por tanto, es movimiento uniforme cuando los parámetros del flujo, concretamente 


la energía por unidad de peso (H), en secciones diferentes, siempre permanecen 
constantes. La forma algebraica de este concepto es: 


H,=H, 


2 2 


v v 
Ad, EOS HL (Ec. 1) 


13 


Diseño de canales abiertos 


En donde: 
H = energía por unidad de peso (m) 
z = energía de posición (m) 
9 = ángulo de inclinación del canal 


a = coeficiente de Coriolis 


Si la altura “H” permanece constante en cualquier sección analizada del cauce, 
entonces en cada tramo, la línea de energía, el plano de la superficie libre y la 
pendiente de la solera del canal, tienen que ser paralelos; siendo la definición de 
movimiento uniforme. Siendo i=I=1, en el movimiento uniforme la fuerza 
predominante para establecer el flujo es la “fuerza componente de la gravedad”, 
actuando paralela con la pendiente del lecho; sin embargo, pueden estar presente 
también las fuerzas netas de presión y las fuerzas de inercia. 


Los flujos uniformes estacionarios ocurren cuando las fuerzas que ocasionan el flujo 
son exactamente balanceadas sobre el tramo en consideración. Este tipo de flujo es 
análogo con el flujo estacionario presurizado en tuberías de diámetro constante. De 
esta manera, el área de flujo uniforme en el canal debe permanecer constante con 
la distancia, y para lo cual es requerido que la pendiente del lecho y la geometría del 
canal permanezcan también es constante. La superficie del líquido es paralela con 
el lecho. 


Finalmente, ensayos experimentales muestran que (% varía entre 1,03 y 1,36 para 
canales prismáticos. Por lo que, 4% >1 por definición (Chow, 1959). En muchos 
casos se justifica Y =1. Sin embargo, es recomendable calcular el valor de alfa. De 
lo contrario siempre habrá incertidumbre en los resultados (% de error). En la Tabla 
1 se resume algunos valores de ( : 


Tabla 1. Valores del coeficiente de Coriolis. 


ES Valor de (Y 
| Mínimo | Promedio | Máximo | 


Canales regulares, canaletas y vertederos 


Corrientes naturales y torrentes 
Ríos bajo cubiertas de hielo 
Valles de ríos, inundados 


1.2. Ecuación del flujo uniforme en canales 
1.2.1. Ecuación de Chézy 


En 1769 el ingeniero francés Antoine Chézy desarrolló probablemente la primera 
ecuación del flujo uniforme en canales. Esta ecuación se puede obtener mediante 
un balance de fuerzas que ocurren en un elemento fluido, no sometido a acciones 
de aceleración. Consideremos la sección de canal mostrada en la Figura 2: 
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Figura 2. Flujo uniforme en canales abiertos. 
Fuente: (Chow, 1959). 


Cuando el movimiento de un fluido es uniforme (v =cte) la aceleración es cero ( 
dv /dt=0) y por consiguiente las fuerzas inerciales también son iguales a cero. 
Por lo tanto, en la “condición de equilibrio” del sistema, las fuerzas actuantes son: 
el peso, la fuerza de presión del fluido y las fuerzas de resistencia al movimiento. Es 


decir: 
Y F =0 
w_-E+f=L =0 
w -f +pA -p,A,=0 

Las secciones transversales analizadas son iguales (A, = A, = A), separadas entre 
sí por una distancia L y perpendiculares a la dirección del flujo. Por tanto, las 
fuerzas hidrostáticas ( £ y £, ) también serán iguales, cuya magnitud es el producto 
de la presión hidrostática ( p= yh ) por el área de la sección (A ). Estas dos fuerzas 


actúan en sentido contrario y a h/3 de la base del canal. En consecuencia, estas 
fuerzas se anulan entre sí. Es decir: 


w -f =0 


La componente efectiva de la fuerza gravitacional (peso) en la dirección del 
movimiento (w_ ) es paralela al fondo del canal e igual a: 


w,=yALsen0 
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En donde: 
w, = componente del peso que origina el movimiento uniforme 
Y = peso específico del agua 
A = área hidráulica de la sec ción del canal 
L= longitud del tramo de análisis 
V = volumen del tramo de análisis 
9 = ángulo de inclinación del canal 


Para valores pequeños de O (generalmente utilizados en canales) sen0 = tan0 , 
por lo que, la componente del peso que origina el movimiento uniforme será: 


w,=7yALtanO 
También, para movimiento uniforme (ver Figura 1): 
¡=/=i =tan0 
En donde: 
¡., = pendiente de la solera del canal 
Por tanto: 


w =yALi, 


La otra fuerza que completa la condición de equilibrio es la fuerza de fricción 
producida entre las paredes y el fondo del canal y el fluido. Así pues, la existencia de 
un gradiente de velocidad implica la existencia de esfuerzos de corte o rozamiento. 


Las fuerzas de resistencia al movimiento corresponden principalmente a la acción 
y reacción de las paredes de la solera del canal y las mismas capas del fluido. Chézy 
supuso que la fuerza que resiste el flujo por unidad de área del lecho de la corriente 
es proporcional al cuadrado de la velocidad (v). Si la superficie de contacto del 
flujo con el lecho de la corriente es igual al producto del perímetro mojado (P ) y la 
longitud del tramo del canal (L), la fuerza total de resiste al flujo ( f) será igual a: 


f e y? 
f =Ky PL 
En donde: 
f, = fuerza de resistencia al movimiento 
K =cons tante de proporcionalidad 
v = velocidad del fluido 
P = perímetro mojado 


L= longitud del tramo de análisis 
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Volviendo a la condición de equilibrio en X : 
YF =0 
Kv*PL-yALi =0 


En flujo uniforme la componente efectiva de la fuerza gravitacional que causa el 
flujo debe ser igual a la fuerza total de resistencia. Es decir: 


Kv*PL=yALi, 


A 
Kv?=y“i 
15 E 
_A 
Pp 
Kv” =yYRi 
YRi, 
V= 
Kk 


En donde: 
v= velocidad del fluido 
Q = caudal que circula por el canal 
A = área hidráulica del canal 
¡, = pendiente de la solera del canal 
C = coeficiente de Chézy 


R= radio hidráulico 
1.2.2. Coeficiente de Chézy 


El número C es un coeficiente que recibió el nombre de coeficiente o factor de 
resistencia de Chézy, R esel radio hidráulico del canal y 1 es la pendiente de la línea 
de energía que para el caso de flujo uniforme estacionario es igual a la pendiente del 
fondo del canal 1 (ver Figura 1). 


La fórmula obtenida por Chézy en 1789 tiene una gran importancia en la historia de 
la hidráulica, ya que de ella se derivan todas las fórmulas modernas para el cálculo 
del flujo uniforme. Sin embargo, quedaba aún el problema de saber cómo valuar la 


17 


Diseño de canales abiertos 


constante C. Chézy solo dijo que ese valor se encuentra entre 30 y 50. 


El coeficiente de Chézy C, está determinado con base a las expresiones matemáticas 
que a continuación se detallan. 


nn” Dr” E , ] Y. a 
1.2.2.1.Ecuaciones analíticas 


El equilibrio que ocurre en un elemento fluido, no sometido a acciones de 
aceleración, considera únicamente dos fuerzas: 1) la componente del peso en 
dirección del movimiento yALsen60, y 2) la fuerza de resistencia al movimiento 
Kv”PL, descrita según Chézy. Es decir: 


% . 
Kv“PL-yALi, =0 
Si el agua fluye en el canal desarrolla una fuerza de fricción en su solera (lecho o fondo) 
y paredes, en dirección del flujo pero sentido contrario. Dicha “fuerza de fricción o 
rozamiento” entre el fluido y el contorno sólido del canal es equivalente a la “fuerza de 


resistencia al movimiento” o “fuerza tractiva unitaria*”. Entonces, la ecuación anterior 
se puede escribir en términos del “esfuerzo promedio de corte (7, ””, es decir: 


2 
T, =Kv 
Resultando: 
T,PL—-yALi =0 
T,PL=yALIi, 
En donde: 


T, = esfuerzo cor tante promedio 

P = perímetro sólido que moja el fluido en la sección transversal 
A= área hidráulica de la sec ción del canal 

R=radio hidráulico en la sección transversal 

Y = peso específico del agua 

p = densidad del agua 

L= longitud del tramo de análisis 


¡ = pendiente de la solera del canal 


Despejando 1, tenemos: 


Como la pendiente de la línea de energía es paralela a la pendiente del fondo del 


1. Llamada también “Fuerza de Arrastre”. 
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canal y a la pendiente de la superficie libre (1= 1 = [), se tiene que la pérdida de 
energía en la distancia L es: 


De donde: 


Pero las pérdidas de energía también se pueden relacionar con el factor de fricción 
(f) de Darcy — Weisback?, cuya expresión es: 


En donde: 
v= velocidad media del flujo 


g = aceleración de la gravedad 


Igualando las ecuaciones (7) y (8) y despejando el esfuerzo cortante (T, ) tenemos: 
T 2 
LÍ_- L yv 
pgR  —4R2g 
RL 2 
To = ¿Pene Y 
4RL 2g 


2 


Igualando las ecuaciones (9) y (4): 


2 
v 


pf--= pgRi, 
8 
Finalmente, la velocidad media en un canal (v) con flujo uniforme será: 
2 8gRi, 
v” = — 


f 


2. Julius Ludwig Weisbach fue un matemático e ingeniero alemán conocido por completar el trabajo de 
Darcy sobre pérdidas de carga en tuberías para dar lugar a la ecuación de Darcy — Weisbach. 
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En donde: 

f =coeficiente de resistencias por superficie 
La dificultad en el uso de esta ecuación radica en determinar el factor de fricción ( 
f ), pues este parámetro depende del número de Reynolds (Re) y la geometría de 
los elementos de rugosidad de los contornos del canal. No obstante, la variable f se 
puede determinar mediante las siguientes expresiones: 


Tabla 2. Factor de fricción de la ecuación Darcy-Weisbach. 


Régimen Fórmula Autor Ecuación 


Pared lisa . = 2log(ReW/f)- 0,8 Prandtl 


Pared de transición- 1 Colebrook y 


21 
régimen turbulento f White 


Pared rugosa Nikuradse 
2k 


Fuente: Elaboración propia. 


En donde: 
Re = número de Re ynolds 
R =radio hidráulico 
k= rugosidad de Nikuradse 


En general, este procedimiento de cálculo se utiliza para cualquier tipo de régimen 
de flujo, pero existe la dificultad en la determinación del coeficiente de rugosidad k . 


Un enfoque sólido del coeficiente C es utilizar la ecuación (12) con un valor de f 
estimado del Gráfico del factor de fricción de Moody. De hecho, para análisis de 
canales abiertos se recomienda usar el factor de fricción f en todos los cálculos. 
Dado que los canales típicos son grandes y ásperos, generalmente se utiliza el flujo 
turbulento totalmente rugoso (White, 1998). 


1.2.2.2. H Xpresiones e 


mpíricas 


Son desarrolladas bajo la consideración que C= £(k, Re, R) . 


Para la mayoría de los casos prácticos, se tiene que el flujo se desarrolla en números 
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de Reynolds lo suficientemente altos (régimen turbulento), en donde el coeficiente 
C deja de depender de éste, por tanto C= £(k, R) 


a) Formula de Ganguillet y Kutter 


En 1869 los investigadores suizos Emile-Oscar Ganguillet y Wilhelm Rudolf Kutter 
obtuvieron una expresión para determinar el valor de C, en función del tipo de 
material y de otras características del flujo. La fórmula de Ganguillet y Kutter es la 
siguiente: 


En donde: 
a=23 
[=1 
m =0,00155 


El valor n es un coeficiente que depende exclusivamente del material y siempre se 
da en el sistema métrico, tal como se determinó originalmente por ser dicho sistema 
el vigente en Suiza. En la Tabla 3 se presenta algunos de los valores más comunes 
del coeficiente n: 


Tabla 3. Valores de n para las fórmulas de Ganguillet - Kutter?. 


Vidrio, tuberías nuevas de estaño, plomo y hierro galvanizado | 0,007- 0,008 


Hierro forjado sin fisuras y tubería nueva de hierro fundido, 


q : Ebo ' . 0,008- 0,009 
recubierto en la mejor condición y alineamiento 

Tubería nueva de hierro fundido, esmaltado y tubería nueva 
acristalada de todo tipo; madera cepillada y en alineamientos | 0,009- 0,010 


perfectos 


Tubería nueva de hierro forjado remachado de diámetro 
pequeño; tubería nueva de duela de madera; madera | 0,010- 0,011 
cepillada, cemento puro 


Madera sin cepillar cuidadosamente unida; cemento, un 
tercio de arena; ladrillo nuevo bien puesto, cuidadosamente 
señalado y raspado; tubería de hierro fundido limpio con 
algún tiempo de uso 


0,011- 0,012 


3. Recuperado de: 
https://archive.org/stream/hydraulicdiagram0Oswanrichittpage/6/mode/2up 
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Madera sin cepillar; cemento, dos tercios de arena; 
mampostería bien cubiertas; tubos de barro y piedra en buenas 
condiciones, pero no es nuevo; yeso y madera cepillada de 
calidad inferior; tubería acristalada mal distribuida o falta de 
uso; nueva tubería de hierro forjado remachado con muchas 
articulaciones y remaches 


0,012- 0,013 


Ladrillo de cara áspera; mampostería bien diseñada y 
distribuida, poco deteriorada por el uso; tubería sucia oO 
ligeramente tuberculada de hierro fundido; tubería grande de 
hierro forjado remachado, pocos años de uso, pero en buen 
estado; revestimiento de lona 


Mampostería en estado inferior o muy sucia; tubería de hierro 
tuberculada y sentida; escombros en cemento o yeso, en buen 
estado; canales de grava rayada con 8,3 pulg., granos bien 


apisonados o rejuntados de cemento 


Escombros en cemento de calidad inferior; escombros gruesos 
fraguados en seco; mampostería en mal estado; canales de 
grava alineada, con granos de una pulgada, bien apisonada o 
cemento con lechada 


Escombros en bruto en malas condiciones; canales con camas 
de tierra en perfecto orden y alineamientos 


Canals con camas de tierra en buen estado y alineamientos y 
libre de piedras y malas hierbas 


Canales con camas de tierra en buen estado de forma 
moderada y alineamientos, con algunas piedras y malas 
hierbas 


Canales con camas de tierra en malas condiciones y 


alineamientos, con piedras y malas hierbas en gran cantidad 
Fuente: (Swan y Horton, 1905). 


La ecuación de Ganguillet y Kutter se dedujo en detalle a partir de datos de 
mediciones de flujo en canales de diferentes tipos, incluyendo los trabajos de Bazin, 
varios ríos europeos y el río Mississipi. La ecuación parece complicada, pero por lo 
general produce resultados satisfactorios. Se ha utilizado con bastante frecuencia y 
existen diversos cuadros y tablas para su aplicación. 
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b) Fórmula de Manning 


En 1889 el irlandés Robert Manning propone la siguiente expresión!: 


En donde C se expresa en (m*?/s) y el número n es conocido como factor o 
coeficiente de rugosidad de Manning. Su valor es relacionado a tipo de superficie 
del contorno. De forma experimental se ha demostrado que el valor n varía con el 
radio hidráulico o la profundidad del flujo. Sin embargo, para fines prácticos se toma 
valores constantes de acuerdo a las características del contorno. En la Tabla 4 se 
listan algunos materiales con sus respectivos valores medios de n: 


Tabla 4. Valores de las rugosidades 1 de Manning?. 


4. F.M. Henderson atribuye a Gauckler (1868) y a Hagen (1881), la observación de que C es proporcion- 
al a R'*, quienes la obtuvieron independientemente; y fue Flamant quien equivocadamente la atribuyó 
a Manning. 

5. Recuperado de: 

http://heidarpour.iut.ac.ir/sites/heidarpour.iut.ac.ir/files//u32/0pen-chow.pdf 
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c) Fórmula de Bazin 


En 1897 el francés Henri-Émile Bazin propone calcular C' con la siguiente fórmula: 


La fórmula de Bazin no es muy usada actualmente para determinar el coeficiente 
de Chézy C. El coeficiente y es un parámetro que depende de la rugosidad de la 
pared, no es el peso específico del agua. Algunos valores se indican en la Tabla 5: 


Tabla 5. Valores de Y para la fórmula de Bazin. 


La ecuación de Bazin se desarrolló primordialmente en canales experimentales 
pequeños. Su aplicación generales menos satisfactoria que la ecuación de Ganguillet 
y Kutter. Existen estudios comparativos entre los coeficientes C' de Chézy, Y de 
Bazin y n de Kutter para los datos experimentales de Bazin y varias corrientes 
naturales. Los resultados demuestran que la ecuación de Bazin no es tan buena 
como la de Kutter, inclusive para sus propios experimentos. 


d) Fórmula de Forchheimer 


En 1923 el ingeniero e hidrólogo austríaco Philipp Forchheimer propone: 


43<C<143, en unidades inglesas. 
e) Fórmula de Pavlovskii 


En 1925 Nikolai Nikolaevich Pavlovskil propone variable el valor exponencial del 
radio hidráulico y determina C de la siguiente manera: 


Para R<l: y=1,54/n Para R>1: y=1,34/n 
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f) Fórmula de Agroskin 
En 1949 el investigador |. Agroskin propone: 


CS as 17,72 log R]|(Ec. 16) 
n 
g) Fórmula de Al'tshul 


En 1952 Adol'f Davidovich Al'tshul' propone para cauces naturales, cuyos lechos 


estén formados con grava y arena, canales no revestidos que acarreen azolves, 
recomiendan la siguiente formula: 


1.2.2.3. Ecuación de Manning 


Reemplazando la ecuación (12) en (2) tenemos: 


Para el caudal tenemos: 


Las Ecuaciones 18 y 19 son las fórmulas que se han convertido en el principal 
instrumento de cálculo práctico para flujos uniformes en canales abiertos. 


25 


Diseño de canales abiertos 


1.3. Secciones hidráulicas 
1.3.1. Sección trapezoidal 


Las variables geométricas en un canal trapezoidal son: el ancho de la base (b); la 
pendiente del talud lateral del canal (m = cot $ ), y el ángulo o inclinación de las 
paredes laterales del canal con respecto a la horizontal (B, Á B,) 


Figura 3. Condiciones geométricas de los canales trapezoidales. 
Fuente: Elaboración propia. 


Área hidráulica: 


b+(m h+b+m,h) % 


Perímetro hidráulico (mojado): 


P=,/h "+(m,h +b+ hi4 (m,h)' 


a qn 
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Radio hidráulico: 


pd, 


— b+h(/L+m; + Lem.) 


R 


Velocidad: 


2 
ld, sl 
1 2 ES 
EA 


.1 2 


Ml b+h(yLem + fem?) j 


Caudal: 


Cuando: m, =m, = 


A=(b+mh)h 


P=b+2hV1+m? 
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(b+mh)h 


Ra == 
b+2hV1+m” 
2 
3| (b-+mh)h l E 
b 


v= 12 
n [b+2hYl+m? | 
5 
as. Pato 
3 


n 
(b+2hV1+m*) 


Para canales triangulares, b=0 (m, 4 m, ), tenemos: 


AE (m, +m,)h 
2 
P=h( Lem] + Lem.) 
(m, +m, )h 


e L (m, +m, )h 
n Alem? + Jm.) 
a | 


l 
.12 
Ss 


2 l 
62 
S 


0=2 


[fama] 


ES 


Cuando: m, =m,=m 
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mh A 


1 
V==:" 
12 


1+m? 
1.3.3. Sección rectangular 


Para canales rectangulares, m, =m,=0 (f, =P, =90"), tenemos: 
A=bh 
P=b+2h 
bh 


1.3.4, Secciones cerradas 


Los canales de sección cerrada comprenden los túneles, conductos y diferentes 
tuberías a través de los cuales el líquido fluye a gravedad. En este caso el líquido 
puede o no llenar toda la sección. Las secciones hidráulicas más usadas son: a) 
circular, b) rectangular, c) de baúl, d) herradura, e) ovalada, entre otras. 


a) 


Figura 4. Secciones de canales cerrados. 
Fuente: Elaboración propia. 
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Para el cálculo hidráulico de canales de sección cerrada se utiliza, también, las 


ecuaciones de Chézy y Manning con algunas adaptaciones, las cuales se detallan a 
continuación. 


1.3.4.1. Conductos circulares de sección llena 


Los parámetros y relaciones hidráulicas más importantes para los conductos 
circulares de sección completamente llena son: 


Área hidráulica: 


D? 
SE NIE 
Perímetro hidráulico: 
P, = 11D 
Radio hidráulico: 
D 
R,=— 
E= Al 
Velocidad: 
0 
Ma DARLA 
¿ln 
43 
Caudal: 
Q,¿=v'A 
Sa 
Or E cid 
43)” 
En donde: 


Aj? área hidráulica del fluido en sección llena (m?) 
PA perímetro hidráulico del fluido en sección llena (m) 
Ry? radio hidráulico del fluido en sección llena (m) 

v y = velocidad del fluido en sección llena (m/s) 

Q, = caudal del fluido en sección llena (m*/s) 


D = diámetro (m) 
¡. = pendiente del canal (m/m) 
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- o O E . , 7 . A A 
1.3.4.2. Conductos circulares de sección parcialmente llena 


Para determinar los parámetros hidráulicos y relaciones hidráulicas de los conductos 
circulares de sección parcialmente llena, introducimos el concepto de ángulo central 
que demarca el sector circular y a la vez representa la zona ocupada por el caudal 


OQ: 


Figura 5. Sección hidráulica parcialmente llena. 
Fuente: Elaboración propia. 
El ángulo central O (en grados sexagesimales): 


0 adyacente 


cos =-— 
2  hipotenusa 
D 
il 
cos-=2— 
2 R 
0 d 
cos = —-=— 
2 2R R 
0 2d 
cos =1-— 
D 
2d 0 
—=1-cos> 
D 
d: 00 
D 2 
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l-cos— 
d= -D 
De lo cual: 
D= 2d 
0 
l-cos— 
2 
Y también: 
2 
6=2 arcos 1 
D 
Si consideramos un franco libre sobre el nivel de agua (S ), tenemos: 
1 0 
D-S cos 
D 2 
EM ad 
D 2 
De donde: 
D= 28 
0 
l1+cos— 
2 
Área hidráulica: 
R? 
A, =(9-seng) 
D? 
Al =(0-sen6).- 
D? 
A, = 0/1 m2) . E 0 = radianes 


Mr 0 -—; O = grados 


TO [- 180 02) D? 
4 


Perímetro hidráulico: 


P, =0-R 
-D ] 
P, = 2D ; O=radianes 
2 
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TO 
P =——D; O = grados 
*=360 a 


Radio hidráulico: 


A 
R, = - 
(9-sen8) —— 
R, = 
9-D 
2 
R, =(9-sen8). 


R, = acne E 0 = radianes 
0 4” 
1 D 
K,=|L= SA ; O = grados 
TO qe 


Velocidad: 
1 LE 
ds 
1 180sen0 | D | > 
V, == 1- 1 
n TO 4 
2 2 
1 180sen0 3 (DY? > 
vi -| 1 . “1 
n TO 4 
2 
3 2.1 
[1 een0] En a 9 = grados 
T n 
43 
2 1 
Caudal: O RP ¡ 
n 


2: 2 
2 3 a A 
Q=% 160 =«end l, DOI ¡ 180sen8 |? D AE 
n 1180 8 TO 4 d 
da 
TO 180sen0 180sen0 Y3 1 ( D? 
180 TO ni 8 


TO 
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5 
3 8 1 
a -[22)( nono $ S —-Di+52: O=grados 


En donde: 
0 = ángulo central (en grados) 
A, = área hidráulica del fluido en sección parcialmente llena (m”) 
1 = perímetro hidráulico del fluido en sección parcialmente llena (m) 


R, = radio hidráulico del fluido en sección parcialmente llena (m) 


v, = velocidad del fluido en sección parcialmente llena (m/s) 
Q, = caudal del fluido en sección parcialmente llena (m?/s) 
D = diámetro (m) 

1 = pendiente del canal (m/m) 


nes hidráulicas en conductos circl 


Relación A, CA : 


TO . 180 me) D? 


A, 360 TÓ 4 
Ay a 
4 
A, _ 0 . 180 0) 
A, 360 TO 
Relación P, /P,, : 10 
P, 360 — 
Pa TD 
P. 08 
P, 360 
RelaciónR, /R,,: 
180sen9 | D 
R, TO ] 4 
Bo Lil 
4 


R,_ 1 180sen0 
Es: TO 
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2 
3 2 1 
(1 ose] E a 
Y TO 43 n 
NA 2 1 
" o 
43 2 
e ¡  180sen0 3 
Ve TO 
Relación Q,/Q, 5 
Oia 1 za 
360 TO ln ; 
Q,_ 43 
Sa E de 
5 Ss 
4)” 
5 
Q...0 ¡ 180 send 3 
Q, 360 TO 
0 E — EA — ' sl |] 


Gráfico 1. Relaciones hidráulicas de una sección circular. 


Fuente: Elaboración propia. 
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Los valores representados en el Gráfico 1 se encuentran en la Tabla 6: 


Tabla 6. Valores de las relaciones hidráulicas — Sección circular. 


R,/R, Y MM, Q,/Q, 
8 


0,000 0,51 | 0,513 | 0,506 
0,000 0,52 | 0,525 | 0,513 7 
0,005 | 0,090 | 0053 | 0,141 | 0,001 0,53 | 0,538 | 0,519 | 1,037 0,551 


0,009 | 0,11 


Y Ss 
p 
05) 


025 016 | 0,534 


a/a, | 0/0 | aa, 049, 


1 
8 
4 
8 
0 
3 
4 


o 
o 
p 
(09) 


0,1 


h 
/0 
0,184 | 0,002 0,54 0,525 | 1,048 0,568 
0,222 | 0,003 0,55 0,532 0,586 
0,257 | 0,005 0,56 0,538 | 1,070 0,603 
0,289 | 0,007 0,57 0,545 | 1,081 0,620 
0,319 | 0,010 0,58 0,551 0,637 

, 0 10 

10 

,0 

,0 


2 
0,14 


0,15 
0,17 
0,18 


0,205 | 0,254 | 0,401 | 0,021 0,61 | 0,639 | 0,571 | 1,078 | 0,689 
0,215 | 0,278 | 0,426 | 0,025 0,62 | 0,651 | 0,577 0,706 


0,225 | 0,302 | 0,450 | 0,031 0,63 | 0,664 | 0,584 | 1,137 0,723 
0,235 | 0,325 | 0,473 | 0,036 0,64 | 0,676 | 0,590 | 1,145 0,740 
0,244 | 0,349 | 0,495 | 0,042 0,65 | 0,688 | 0,597 | 1,1 | 1,099 | 0,756 


0 

0,194 | 0,230 | 0,375 | 0,017 0,626 | 0,564 111 | 1,072 | 0,672 
h 7 

0,253 | 0,372 | 0,517 | 0,049 | 0,700 | 0,604 104 | 0,773 

11 

12 


0,205 | 0,348 | 0,013 0,59 0,558 | 1,101 0,655 


0, 
0,538 | 0,056 0,67 0,610 | 1,108 | 0,789 
0,558 | 0,063 0,68 0,617 | 1,173 0,806 
0,279 0,577 | 0,071 0,624 | 1,179 0,821 
0,287 0,597 | 0,079 0,70 0,631 | 1,185 | 1,120 | 
0,295 0,615 | 0,088 0,71 0,638 ,12 
0,303 0,633 | 0,097 0,72 0,645 
0,311 0,651 | 0,106 0,73 0,652 | 1,199 0,883 
0,318 0,668 | 0,116 0,74 0,659 | 1,203 0,898 
0,326 0,684 | 0,126 0,75 0,667 | 1,207 0,912 
0,333 0,701 | 0,137 0,76 0,674 
0,717 | 0,148 0,77 0,682 y 
0,732 | 0,159 0,78 0,689 0,953 


o 
o 
00 
¡05 
N 


s ES 

o 

u 
El 
o 
00 
a 
w 


un 


z 


158) 
m 
o 


7 


o 
un 
o 
lo] 
155] 
m 


o - 5 - 
e) 00 
Y D 
y 00 


IS) 
05] 
N 
o 
lo] 
uy 
[Co] 


m 


3 
3 
0,28 0,646 | 0,747 | 0,171 0,79 | 0,847 | 0,697 0,965 
0,29 0,665 | 0,762 | 0,183 0,80 | 0,858 | 0,705 4 


0,30 0,684 | 0,776 | 0,196 0,81 | 0,868 | 0,713 


0,31 0,790 | 0,209 0,82 0,721 1,000 
0,32 0,804 | 0,222 0,83 0,729 

0,817 | 0,235 0,84 0,738 1 
0,830 | 0,249 0,85 0,747 187. | 9, 
0,774 | 0,843 | 0,263 0,86 0,756 1,039 


(NS A 
o 


un 


o 
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0,39 | 0,361 


Y, 


d/D 
0,87 
0,88 
0,89 


0, 


0,924 | 0,765 


0,932 


0,940 | 0,785 


0,948 


0,775 


0,795 
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1,207 [1 


H h H 

0,410 

0,416 

0,423 

; 0,429 
0,449 
0,455 
0,462 
0,468 
0,475 0,97 

0,98 


1,00 | 1,000 


0,806 
0,817 


0 
Posa | 0:13 | 0455 | 0505 | 05m 
Pass | 0430 | 0462 | 0558 | 03m 
5 


0,93 0,830 


1 0,94 0,842 


0,41 
0,424 0,95 0,856 1,095 
2 | 0,48 

4 


0,872 1,086 


0,889 1,075 


0,910 
0,936 


0,494 1,000 


No obstante, debido a la gran complejidad de las expresiones algebraicas para la 
determinación del área, el perímetro mojado y el radio hidráulico, también suelen 
utilizarse Tablas o monogramas para simplificar el cálculo. Entre los métodos de 
mayor sencillez tenemos los Gráficos o ábacos, en el cual se utilizan las relaciones: 


a=f| —|f Q, 

K, Qu 
E 
Vy 


En dónde: 

a = relación de llenado para el caudal dado. 

K,, Q,= factor de gasto y el caudal para el tirante d . 
K., Ds = factor de gasto y el caudal para la sección llena. 
v, = velocidad del flujo con el tirante d . 


eS velocidad del flujo para la sección llena. 
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En el siguiente gráfico se presenta las relaciones hidráulicas para conductos de 
sección tipo baúl. 


Q, =2,2164(/1/n) (5/2)*” ; A, =3,544(H/2)* 


h/H 


TIPO BAUL VH/ Vi 


Gráfico 2. Curvas de caudal y velocidad de conductos tipo baúl. 
Fuente: Elaboración propia. 
1.3.5, Resumen 
Afin de facilitar los cálculos hidráulicos se resumen las propiedades de los diferentes 
tipos de canales abiertos. 
Para canales trapezoidales, triangulares y rectangulares: 


Tabla 7. Ecuaciones para canales abiertos. 


Parámetro | Fórmula | Ecuación 


Cuando 1, 4 'M) (canal trapezoidal) 


Área 
hidráulica 


Perímetro 
hidráulico 
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Radio 
hidráulico 


Velocidad 


y= 


O= 


Cuando M, = M, = M (canal trapezoidal) 


An 


_— bh fem? +fLem]) 


n ph fem? + f1+m,) " 


> [pel lenin) 


1 
. [2 
2% 
3 


Pablo Gallardo Armijos 


[? 


Área Ñ 
o 


Perímetro 
hidráulico 


Radio 
hidráulico 


Velocidad 


P=b+2hWV1+m? 


(D+mh)h 


R = 
b+2hV1+m? 


(D+mh)h 


1 e 
=p AA á 
[(0+mn)a] 
AAA 


1 
g= 7 Í, 
[6+20V1+m) 
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Área 
hidráulica 


Perímetro 
hidráulico 


Radio 
hidráulico 


Velocidad 


Cuando b= 0 y M, 4 M, (canal triangular) 


P=hf,[+m; + tem) 
(m, +m,)h 


2 Lem: +ftem?) 


R= 


(m, +m,)h 


2 Lem: + fLem]) 


(m, +m,) 1 ; 


O= E 
al 


E 
v==—+ pe 
n 


1 
.p2 
e 
3 


Área 
hidráulica 


Perímetro 
hidráulico 


Radio 
hidráulico 


Velocidad 


Cuando b=0 y M, = MM, (canal triangular) 


A=mh? 


P=2h/1+m? 
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(mie) 


e 
;52 


021 
3 


Ñ (24m? 


Cuando b=0 y M,= M, (canal rectangular) 


Área 
hidráulica 


A =bh 


Perímetro 
hidráulico 


P=b+2h 


Radio 
hidráulico 


Velocidad 


Caudal 


Fuente: Elaboración propia. 
Para conductos de sección circular: 


Tabla 8. Ecuaciones para conductos circulares. 


Parámetro | Fórmula | Ecuación 


Parcialmente llena 


á 2d 
AnEnIa 0=2 arccos| 1-77) (Ec. 45) 
central D 


0 
Tirante _ 2 » 
hidráulico ¡sii 
2 
Área a 91 Meno) 1 
hidráulica 360 4 
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Perímetro 
hidráulico 


Radio _f, 180seng] D 
hidráulico TO 4 


2 
| Se) 
Velocidad LA 1— : 
TO 

5 

| 

Caudal = | 1-——— |. 
TO 


Área 
hidráulica 


Perímetro 
hidráulico 


Radio 
hidráulico 


Velocidad 


Caudal 


Fuente: Elaboración propia. 
Las relaciones hidráulicas en conductos circulares son: 


Tabla 9. Relaciones hidráulicas de conductos circulares. 


Áreas 1 180 2) 


hidráulicas TO 
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Perímetros 
mojados 


Radios 
hidráulicos 


Velocidades MN 21. a 


5 
Caudales 24 :z 0 pp | 


O, 360 TO 


Fuente: Elaboración propia. 


1.4. Determinación del calado normal 


Se define como caldo o tirante normal a la profundidad que se establece en un canal 
con flujo uniforme. 


Existen varios métodos para el cálculo del calado. Los más conocidos son dos: 
+ Método del factor de gasto o caudal característico 
+ Método del exponente hidráulico 


Estos métodos son aplicables para cauces naturales y prismáticos. Conociéndose 
como cauces prismáticos aquellos que tienen sección transversal de forma 
geométrica definida (trapezoidal, rectangular, triangular, circular, etc.). 


Para el análisis respectivo son necesarias las variables geométricas descritas en la 
Figura 3, es decir: ancho de la base (b ); pendiente del talud lateral (m ), ángulo de 
las paredes laterales ( 8 ), caudal (Q ) y pendiente de la solera (1, 1: 


1.4.1. Método del factor de gasto o caudal característico 


Q=CA[Ri, 
(Ec. 61) 


K = factor de gasto o caudal? (m?/s) 


De la ecuación (3): 


Suponemos: 


En donde: 


C= coeficiente de rugosidad de Chézy 


6. Llamado también conductividad de una sección de canal. 
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A = área hidráulica del canal (m?) 
R = radio hidráulico del canal (m) 


Según la ecuación (61) el parámetro K depende directamente de la rugosidad y de 
las condiciones geométricas del canal. No obstante, al reemplazar la ecuación (61) 
en la ecuación (3), el valor de K también es igual a: 


a 


Ea 

1 
De lo cual, para un caudal (Q ) y una gradiente hidráulica (1,) conocidas, existirá un 
solo valor de K.. En este caso, el factor de gasto o caudal ( K ), se denomina “gasto 
o caudal característico (K, )”, de tal manera que: 


K, = Q (Ec. 62) 
di 
En donde: 


K, = gasto o caudal característico (m?/s) 
Q = caudal de diseño (m?/s) 
ja pendiente de la solera 


El parámetro K, depende directamente de las condiciones hidráulicas del flujo. Por 
lo tanto, dadas las condiciones de caudal (Q ) y gradiente hidráulica (1,) se puede 
conocer directamente el caudal característico ( K, ). 


Finalmente, para determinar el calado normal del canal (h), la solución consiste en 
determinar un valor de (K) de tal manera que K, =K.. Es decir: 
K, =K 


Q_ 
E CAR 


Luego, reemplazando la ecuación (12) EOS 


AR 
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En donde: 
K, = gasto o caudal característico de diseño (m?/s) 
n = coeficiente de rugosidad de Manning 
A = área hidráulica (m?) 
P = perímetro hidráulico (m) 


La ecuación (63) es la fórmula para el cálculo práctico de flujos uniformes en canales 
abiertos. El término de la derecha se denomina “factor de sección (s)” y depende 
únicamente de la geometría mojada del canal. 


Se calculará el factor de sección de ) hasta obtener la siguiente igualdad: 
K, .n=F 


(Ec. 65) 


Esta condición se cumple dándole valores a (h ) en las ecuaciones (20) y (21) (aplicado 
también para secciones circulares) y calculando el factor de sección ( F) hasta obtener 
el valor buscado (K, -1n), por aproximaciones sucesivas. 


Según la forma del canal (trapecio, rectángulo, triángulo o círculo), el factor de 
sección (E) crecerá conforme lo hace el calado normal (h.). En general, existirá un 
único valor de calado normal para que el factor de sección E sea igual al producto 
E. 


Para secciones trapezoidales: 
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Para secciones circulares parcialmente llenas: 


3 
¡180 me. 


TO 


1.4.2. Método del exponente hidráulico 


Según Bakhmeteff los factores de gasto con sus profundidades, para un mismo 
canal, guardan la siguiente relación: 


En donde: 
X = exponente hidráulico del cauce 
K, y K, = factores de gasto o caudal (m3/s) 
h, y h,= profundidades correspondientes a K,yK, 
El exponente hidráulico es considerado constante para cada cauce. 


La solución, según este método, consiste entomar dos profundidades y y determinar 
y, de tal manera que: 


expresión: 


1.5. Canales de contorno rígido 
1.5.1. Sección hidráulica óptima 


Según la ecuación (61) el “factor de gasto o conductividad de una sección de canal” 
se incrementa con el aumento del radio hidráulico o la disminución en el perímetro 
mojado (Chow, 1959), es decir: 


K=CAWVR = ca £ 
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Por tanto, para un área determinada de la sección de canal (A ) que tenga el menor 
perímetro mojado (P ), se tendrá la máxima conductividad (K ). Tal sección se 
conoce como “sección hidráulica óptima”. En otras palabras, la sección hidráulica 
óptima es aquella sección transversal (A) de un canal que tiene la máxima 
capacidad, escurrimiento o caudal. 


Así pues, dentro de todas las secciones conocidas, el semicírculo es la sección que 
tiene el menor perímetro mojado para un área determinada, por lo tanto, es la 
“sección hidráulica más eficiente de todas las secciones” (Chow, 1959). 


Analizando la ecuación (63) notamos que, tanto el área A , el gradiente 1 y la 
rugosidad M de canal, están dados (son datos), por lo que, un valor máximo del 
caudal Q obtendremos para 1: que corresponde a un valor mínimo del 
perímetro mojado Pa . Es decir: 


La solución consiste en buscar el mínimo perímetro mojado para un área hidráulica 
dada. En la Figura (3) tenemos: 


da Ea h 


Y: 
P=b+(/L+m; +,/Lem.)h 
De donde: 
b=2_ m,+m, h 
h 2 


Por consiguiente: 


A (m>+m 
P= Amin), +n( Lem fem) 
A m,+m 
P=7+h Jim rfi 0) 


Considerando una sección A dada y un valor m constante, podemos derivar para 
h y obtener: 


+ 


47 


Diseño de canales abiertos 


Reemplazando A tenemos: 


m +m 
me + firm a fimo az Jo 
A A 2 


b=|,Ltm/ +. 14m, —(m, +m) |h 


Si tomamos la segunda derivada de esta ecuación obtenemos un valor positivo, lo 
que dicha ecuación es el mínimo de la función. 


Con la relación hidráulica óptima de b en función de h (para un canal trapecial), 
podemos obtener el área, perímetro y radio hidráulico óptimos son: 


a=([ tem; +/[L+m,? -(m,+m )po[ 2% Jo 


+ 
a |[yiem +, [Lem en peo ln) 
2 
m,+m 
A= Vin +, /1+m* (Bm ho 
Reemplazando b en la ecuación del haa tenemos: 


P= Jim; + +m,? - (m, +m, )Jn+(JL+m; +. 14m, | 


P=2(,/1+m, +/Lem)h-(m, +m,)h 


p=a] [Jim fimo) (Pon) h 


El radio hidráulico es: 
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+ 
rem +, /1+m, mn) lv 


R = 
A n m,+m 
NEO +.[L+m, mm) 
h 
R=-> 
2 
La Ec .60 quedaría de la siguiente forma: 
AS? 
K, =p 


5/3 
m +m 
l l+m?+,/1+m,? (mum)! 
+ 
El im afirmo (M0) 


+ 
| E] pi 
1 
a 2/3 
+ 
rem from 1050] pos 


De donde: 
3/8 


E 2.K,-n 
m +m 
Jl+m?+,/1+m,* mms) 


En la Tabla 10 se presenta un resumen de los parámetros hidráulicos para canales 
trapezoidales y rectangulares con sección óptima: 
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Tabla 10. Ecuaciones para h y D en secciones óptimas de canales. 


Parámetro | Fórmula | Ecuación 


Cuando ato 


Tirante 
hidráulico 


Tirante 
hidráulico 


Tirante 
hidráulico 


Fuente: Elaboración propia. 


Para determinar los parámetros hidráulicos: A, P, R, v y O se utiliza las 
ecuaciones de la Tabla 7. 


El concepto de la sección hidráulica óptima se aplica solo para el diseño de canales 
de contorno rígido (no erosionables). 


| 


1.5.2, Sección económica 


Desde el punto de vista práctico, la sección hidráulica óptima es la sección que da 
el área mínima para un caudal determinado, pero no necesariamente es la mínima 
excavación (sección más económica). 


La sección económica (sección con mínima excavación) sólo ocurre cuando el nivel 
del agua de la sección óptima coincide con el nivel del terreno. No obstante el nivel 
de agua casi siempre se encuentra debajo del nivel del terreno dejando un “franco 
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de seguridad” o “borde libre””. Todo canal más angosto que la sección óptima tendrá 
la mínima excavación (Chow, 1959). 
El área de excavación mínima será: 
m,+m 
A =|B+| 2 


e 


H|H 


Siendo: 
H=S+h+e 
B=b+ 2, AL, 


El tirante hidráulico h y la base del canal b corresponden a la sección óptima del 
canal, variables que son calculadas con las ecuaciones de la Tabla 10. Por semejanza 
de triángulos tenemos: 


x,=mye 
x,=m,e 
De donde: 
B=b+(m, +m,)Je 
m+m 
A,= b+(m, +m,)e+ 2 H 
2 


B 


Figura 6. Cálculo del área de excavación de un canal. 
Fuente: Elaboración propia. 


7. El borde libre es la distancia vertical desde la parte superior del canal hasta la superficie del agua en 
la condición de diseño. Esta distancia debe ser lo suficientemente grande como para prevenir que las 
ondas o fluctuaciones en la superficie del agua causen reboses por encima de los lados. 
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El volumen de excavación por metro lineal (en m*/m) será: 
V=A Im 
e e 


(Ec. 76) 


Para m, =m,=m tenemos: 


oro) 


Por su parte, el área mínima de hormigón requerida para el revestimiento del canal 
se puede obtener mediante la siguiente expresión: 


A, = (lem? + /1+m,)(5+1)+3(m,+m,) e 


b x 


B 


Figura 7. Cálculo del área de hormigón de un canal. 
Fuente: Elaboración propia. 


De igual forma, el volumen de hormigón será: 
V, =A,lm 


v,= (yema +, [Lem )[S+h)+3(m, +m,)+b [e (Ec. 77) 


Para 1,1, =m tenemos: 


v, = | L+m* )(S+h)+6m+b le 


52 


Pablo Gallardo Armijos 


1.6. Canales de contorno erosionable 


Para el diseño de canales erosionables son comúnmente usados dos criterios 
fundamentales: la teoría de la “fuerza tractiva crítica” y el de la “velocidad máxima 
permisible”. 


1.6.1. Teoría de la fuerza tractiva crítica 
La fuerza ejercida por el agua sobre el área mojada de un canal es llamada “fuerza 


tractiva. De la deducción de la ecuación (4) se puede concluir que: “el promedio de 
la fuerza tractiva unitaria es el promedio del esfuerzo cortante”, dado por: 


T, = PERi, 
En donde: 
Z, = esfuerzo cortante del entorno del canal, N/m? 
P = densidad del fluido, kg/m? 
g = aceleración de la gravedad, m/s? 
R = radio hidráulico, m 
j = pendiente de la solera del canal, decimal 


Sin embargo, el esfuerzo cortante del contorno T, no es uniformemente distribuido. 
La distribución varía un poco más con la forma del canal pero no con el tamaño. 


wh-i, 


Figura 8. Distribución del esfuerzo cortante sobre el contorno. 
Fuente: Elaboración propia. 


De lo anterior, el esfuerzo cortante en la solera del canal (7,) será: 


T, = PERi, 
El esfuerzo cortante máximo en los taludes laterales del canal (7, ) será igual a: 
T,=0,767, 


T_=0,76pgRi | (Ec. 78) 
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El esfuerzo cortante en la solera y los taludes laterales del canal debe ser mantenido 
de bajo de aquel valor crítico (T_) que causará que el material del canal se mueva. 
Entonces el canal será estable. La fuerza tractiva crítica (7,) de una partícula del 
material es la fuerza tractiva unitaria que no causará erosión del material sobre una 
superficie horizontal. 


El material sobre los taludes del canal es sujeto, además del esfuerzo cortante debido 
al flujo del agua, también a la fuerza de la gravedad debido a la inclinación de los 
lados. Se ha demostrado que si Ty ES la fuerza tractiva crítica, el máximo esfuerzo 
cortante crítico debido al flujo del agua sobre los lados del canal es: 


En donde: 
Ty 7 fuerza tractiva crítica sobre el lecho del canal, N/m2 
1.5 fuerza tractiva crítica sobre los laterales del canal, N/m2 
0 = ángulo de fricción o ángulo de reposo del material 
Y = ángulo o inclinación del talud del canal con relación a la horizontal 


Igualando las ecuaciones (78) y (79), y despejando R tenemos: 


El primer valor de h puede ser asumido mediante la suposición que R = h, es decir: 


1 sent O 
_ $ sen*0 
0,76pgi, 


Con el valor de h y mediante un proceso de prueba y error podemos obtener el 
valor de b utilizando la ecuación (63): 


A 
Ki 
p3 
Siendo: 
xk - 
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Quedando la ecuación: 


Q (bh+m?P 


>| (Ec. 81) 


E 
vi, [b+2hV1+m J 


Finalmente, se calcula el valor de R mediante la expresión: 


2 
O bh+mh 


o b+2hV1+m? 


1.6.2. Velocidad máxima permisible 
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La velocidad máxima permisible o velocidad no erosionante es la mayor velocidad 
promedio que no causará erosión en el cuerpo del canal. Esta velocidad es muy 
incierta y variable, y sólo puede estimarse con base en experiencia y criterio. En 
general, los canales viejos que han soportado muchos periodos hidrológicos 
permiten velocidades mucho más altas que los canales nuevos, debido a que su lecho 
a menudo se encuentra mejor estabilizado (Chow, 1959). En la Tabla 11 se resume 
las velocidades máximas permisibles y los esfuerzos cortantes para diferentes tipos 


de suelo: 


Tabla 11. Velocidades máximas y fuerzas tractivas para algunos materiales. 


Tierra negra arenosa 


Tierra negra limosa 


Velocidad | Fuerza | Coeficiente de 
Tamaño media tractiva | rugosidad de 
(mm) aprox. crítica Manning 
(m/s) (N/m2) (n) 
L 


0,020 
0,020 


Arcilla (muy coloidal 
0,062-0,2 


-0, 5 
Arena gruesa (1) 0,50-2,00 


1) [E 

imosa (1) E] 

Limo aluvial (1) || 

aa 
Tierra negra ordinaria 

ani [lA 

idal) ¡IA 

1) | 0,062-0,25 | 


Arena media (1 


0,25-0,50 


1,15 12,2 0,025 
0,75 0,020 


3, 
ñ 
1, 
2, 
3, 


0 
5 
5 
7 
7 
2 
7 
5 
7 


0,020 


[37 | 


0,020 


0,020 
0,020 


| 12 | 
| 17 | 0020 | 
| 25 | 
| 37 | 


Grava gruesa 8-64 1,25 14,7 0,025 
Canto rodado y ripio 64-256 1,55 44,0 0,035 


Canto rodado y tierra negra 


0,004-64 1,15 
graduadas (1) j ! 
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(1) No coloidales, (2) Coloidales. 
Fuente: Elaboración propia. 


A partir del criterio de la velocidad máxima permisible, el procedimiento de cálculo 
es el siguiente: 
e Con base del tipo de material del cuerpo del canal, se determina la 
velocidad máxima permisible (V y el coeficiente de rugosidad (n) 
(ver Tabla 11). 


e Con la pendiente longitudinal (1,) se calcula el radio hidráulico (R 
), mediante Manning (Ecuación 18): 


max ) 


e Conel caudal (Q) y la velocidad permisible ( V ix ) Se calcula el área 
hidráulica (A ), es decir: Q 


e Con el valor del radio hidráulico (R ) calculado en el numeral 2 se 
obtiene el perímetro mojado ( P ) de la relación: 


p=2 
R 


e Con el valor del radio hidráulico (R ), la pendiente longitudinal del 
canal (1,) y el ángulo de reposo del suelo (0 ) se calcula el ángulo 
lateral del talud (YH ). Para ello igualamos las Ecuaciones 78 y 79: 


sent O 

0,76pgRi =7T ,|l— 
PERS = Ta sen?0 
0,76pgRi, Í a sen" Y 
Ta sen?0 


(Ec. 82) 


2 
Ri 
senO =sen60, |1— ES 


En 


e Conocido el ángulo lateral del talud (O ) se puede calcular el talud 
lateral del canal (m ) mediante la igualdad: 


A, 
m 
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De donde: 
m = cotan Y 


e Conlos valoresde A, P y m se obtiene las dimensiones del canal, 
planteando un sistema de dos ecuaciones, ambas en función del 
ancho (b) y del tirante hidráulico (h). El tirante hidráulico (h) 
resulta de la solución de una ecuación cuadrática: 


P = b+h|,/I+m +, 14m 


b=P- h(,JL+m? + 14m, 


a 
a o O 


(m, +m,) +m,) 


(Lem +./1+m,?)h Jr +Ph-A=0 


[25 (rem fi5m) +Ph-A=0 


Kh?+Ph-A=0 
En donde: 


K == (m, +m,) (+ m¿ +. 14m, ) 


La solución de la ecuación tiene la siguiente forma: 


_ —b+yb?*-4ac 


2a 


Entonces: 


—P+/P?+4KA 


2K 


h = 


No obstante, dadas las características geométricas del canal (m, b y h) y 


las condiciones mecánicas del suelo (YH), es posible obtener la velocidad (v) 
correspondiente, mediante la ecuación: 
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ui] 


sen*0 

1 To, l= 20 1 

ya =. WE [12 (56,83) 
n 0,76 pg, ] 


Esta velocidad no deberá superar la velocidad máxima de la Tabla 11. 
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CAPÍTULO 2. FLUJO NO UNIFORME EN CANALES 


2.1. Condiciones generales 


Una de las principales tareas de estudio del flujo no uniforme en canales es la 
determinación de las diferentes formas de los perfiles de la superficie libre y por 
consiguiente los tirantes que se establecen. 


Conocemos que para el flujo no uniforme la velocidad media (v), el tirante (h), la 
pendiente de la superficie libre (h ), etc., no se mantienen constantes a lo largo del 
flujo, Figura 9. 


Para cauces con 1 constante, el tirante h cambia a lo largo del cauce y el perfil de la 
superficie libre es una línea curva. 


Figura 9. Flujo no uniforme en canales. 
Fuente: Elaboración propia. 


Si el flujo es retardo, Figura 10a, el tirante h aumenta a lo largo de éste y la superficie 
libre forma una curva llamada de remanso o presión. Para un flujo acelerado el 
tirante disminuye y la superficie forma una curva de derrame o depresión, Figura 
10b. 


a) Remanso b) Derrame 


Figura 10. Formas de la superficie libre. 
Fuente: Elaboración propia. 
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2.2. Ecuación diferencial del flujo no uniforme 


Recordemos que un flujo gradualmente variado implica que la distribución de la 
presión es hidrostática para el plano normal, al del flujo, lo que nos permite aplicar 
la ecuación de Bernoulli. En la Figura 9 se puede ver que: 


H=z+h 


Diferenciando y dividiendo para —dl obtenemos: 


dH_ dz dh 


Como: 


Resulta que: 
. dh 
I=1.—— (Ec. 84) 
* dl 


Por otro lado, la energía específica entre dos secciones está dada por la ecuación de 


Bernoulli. a 


v 
z+h+0—+h, =cte 
28 
Reemplazando z+h= H, tenemos: 


H+0—+h = cte 
la * 


Derivando para dl : 


De donde: 


Pudiendo escribir también: 
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Además: 


2 
de) 
284 ")_0Q07 14) (Ec. 86) 


dl 2g di 


Como el área de una sección invariable en el tiempo es función de la longitud (1) y 
la profundidad (h), el diferencial total de la función área (A) es: 


JA JA 


dA = —dl+-— dh 
ol oh 
Por lo que: 
1 dA 2 [JA JA 
dl + |=-2—=-—| —dl+ —dh 
(2) A? e dh j] 


(Ec. 87) 


a | J 224,24 4) 


d Aélol ah md 


Un incremento de área JA =0B-h, donde B es el ancho del cual en la superficie 
libre, tenemos que: 


0 pics) 
dh 


Por lo tanto, reemplazando las ecuaciones 87 y 88 en 86 tenemos: 


d av 
28 aQ 2(0A _dh 
AAA — + B— 
dl 2g A“Lol dl 
gh 
dl 
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Igualando con la ecuación 84 tenemos: 
dh 00Q? JA +B dh Q? 
* dl gA? ol dl A?CR 
dh 00? JA aq? SS dh , Q? 
* dl gA? dl ga? dl A?CR 
Ez 0Q7 JA. Q% _dh 00Q gh 
5 ga ol ACR dl ga? dl 
2 2 2 
00 JA Q 2-2 a] 


1 


PS mE 


a 0aQ IA Q? 
dh * | gA? 01) A“CR 


dl Ñ 000B 
ga” 
Esta es la ecuación diferencial del movimiento no uniforme para canales abiertos 


de sección variable. En causes de sección constante, el área Á no depende de la 
longitud 1, por consiguiente dA / dl= 0, tomando la ecuación anterior la siguiente 


forma: Q? 
dh ACR 
dl ia 0.0B 
ga” 


Siendo esta la ecuación diferencial de la superficie libre para canales prismáticos. 
2.3. Energía específica de una sección 


Como energía específica de una sección se conoce a la parte de la energía del flujo, 
determinada por el tirante h y la altura de velocidad v”/2g, sin la consideración 
de la energía específica de posición, de tal manera que: 


Z 
she (Ec. 92) 
28 


62 


Pablo Gallardo Armijos 


Figura 11. Energía específica de una sección. 
Fuente: Elaboración propia. 


Para el flujo uniforme la energía específica de una sección se mantiene constante, 
osea: 


28 28 


Para un canal, con caudal dado, la energía específica de la sección es función del 
tirante h y siempre es positiva (e > 0): 


e=f(h) 


2 
Rs 
2gA (Ec. 93) 


Si analizamos gráficamente la ecuación 93 tenemos que, 


Si h > 0, entonces A > 0. Por tanto : po Soye>o 
g 


2 
l 10 
Si h > oo, entonces A > oo. Por tanto : 3 a >0ye>0 
8 
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La energía específica de la sección se encuentra representada en la Figura 12: 


a... e 
min 


Figura 12. Variación de la energía específica. 
Fuente: Elaboración propia. 


2.3.1. Tirante crítico 


La ecuación 3.96 tiene un valor mínimo, correspondiente a un valor de la profundidad 
llamada tirante crítico (h.). El tirante crítico corresponde a la energía mínima de 


cr 


una sección, y analíticamente se obtiene derivando la Ecuación 92: 


ho 
28 

2 

d+ 

de 2gA 
dh dh 


Para un cauce de sección rectangular es fácil determinar la magnitud del tirante 
crítico, ya que: 


A=Bh 
2 
dl 1+ 22 
2gBh 
0O= 
dh 
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En donde: 


Introduciendo el concepto de caudal unitario o específico a la relación tenemos: 


_Q 
1=B 


El tirante crítico tiene gran importancia física, ya que permite dividir a las corrientes 
en dos tipos de flujo: 


Si el flujo posee una profundidad mayor que el tirante crítico éste se denomina flujo 
subcrítico, tranquilo o fluvial. 


Si su profundidad es menor al tirante crítico se denomina flujo supercrítico, rápido 
o torrencial. 


Es decir: 
h> h,, > subcrítico 


h < h,, > supercrítico 


2.3.2. Pendiente crítica 


Se llama “pendiente crítica” a la cual, dado un caudal Q, se establece un flujo uniforme 
con una profundidad igual al tirante crítico. La pendiente crítica se encuentra a partir 
de la ecuación de Chézy. 
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Q=CAVRi 
Q? 
CAR 
2 
1, = 2 E 2 
Co, A, R. 


Donde Ces el coeficiente de rugosidad de Chezy. Según Manning puede ser calculado 


con la siguiente relación: Rv 
pa cr 
E. > 
n 
Por tanto, para canales rectangulares: 
Q? 
cr 1/6 2 
cr 2 
A Eu 
n 
; no? 
ij R EN 2 
cr cr 
n?P 0 
1 = cr 
cr A 4BA 2 
cr cr 
4/3 
n”(B+2h,) 0” 
=> er 
er 4/3 2 
(Bh,,) (Bh,) 
cr cr 


La pendiente crítica sirve también como parámetro de comparación para la 

determinación de flujos subcríticos y supercríticos en canales. Si la pendiente del 

canal es menor que la crítica el régimen es fluvial y si es mayor es torrencial. Es decir: 
¡<1,, > subcerítico 


1>1.. > Supercrítico 
cr 
2.4. Superficies libres 
2.4.1. Formas de la superficie libre 


Para un valor de pendiente de solera del canal ñ y de pendiente crítica 1. se puede 
decir: 
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1 <1. —> pendiente suave 
Ss cr 
1 >1. —> pendiente pronunciada 
S cr 
1 =1. —> pendiente crítica 
Ss cr 
i¡ =0> pendiente horizontal 
¡ <0 > pendiente adversa 


Las diferentes formas de la superficie libre se pueden determinar según la Ec. 91. Es 


decir: . Q? 
dh ACR 
dl e aQB 
ga? 


67 


Diseño de canales abiertos 


Endonde: K, esel caudal característico del flujo uniforme correspondiente al tirante 
hO0; K es el Factor de gasto de la profundidad existente h; el término aQ /g es 
la función que corresponde a la profundidad crítica (h. ,); y el término A*/b esla 
función que corresponde a la profundidad existente (h). 


La Ec. 97 expresa la variación de la profundidad a lo largo del flujo. Un valor positivo 
de Ah/Al indica curvas de remanso y un valor negativo indica curvas de derrame. 


2.4.2, Diseño de la superficie libre 
Para el análisis de la superficie libre simplificaremos la Ec. 97 a la siguiente forma: 
h 
1-f| 
AB; h 


Al ” 
ES 
h 


La interacción entre tirante normal y tirante crítico, originan en el flujo las formas de 
superficie que se exponen en la Tabla 12. 
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Tabla 12. Formas de la superficie libre. 


Formas de los perfiles Análisis de la Ec. (97) Ejemplos 


Tipo M, 
h>h,>h,, 
F(n,/h)<1 
Fr, /h)<1 
Pendiente suave Ah/Al=+/+ 
SO Remanso 
Tipo M (Mild) Tipo M, 
h,>h>h, 


F(n,/h)>1 


F(h,/h)<1 
Ah/Al=-/+ 


Derrame 


Tipo M, 
h,>h, >h 
F(n,/h)>1 
F(»,,/h)>1 

Ah/Al=-=/- 


Remanso 
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Formas de los perfiles Análisis de la Ec. (97) Ejemplos 


Tipo S, 
h> her > h, 
F(h,/h)<1 
Fr, /h)<1 
Ah/Al=+/+ 


Remanso 


Pendiente pronunciada 
i>i 
Ss cr 
Tipo $ (Steep) 


Tipo S, 
her >h> h, 


F(h,/h)<1 

F(»,,/h)>1 

Ah/ Al=+/- 
Derrame 


Tipo S, 
her > h, >h 


F(n,/h)>1 
F(»,,/h)>1 
Ah/Al=-/- 


Remanso 


Pendiente crítica 


Ah/Al=+/+ 


Remanso 


Remanso 
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Formas de los perfiles Análisis de la Ec. (97) Ejemplos 


Tipo H, 
h>h 


Pendiente horizontal F(n, / h) >1 
¿=0 F(h,/h)<1 
Ah/Al=-/+ 


Derrame 


Tipo H 


Tipo A, 
h.>h 
F(n,/h)>1 
F(»,,/h)>1 
Ah/Al=+/+ 


Remanso 


Tipo Á, 
h>h. 
Pendiente adversa F(n / h) >1 
0 
¡=0 Fr, /h)<1 
Tipo Á Ah/Al==/+ 
Derrame 
Tipo A 
h., > h 
F(n,/h)>1 
Fa, /1)>1 
Ah/Al=-/- 


Remanso 
Fuente: Elaboración propia. 


La Ec. 97 es la expresión diferencial de la forma de la superficie libre y se resuelve 
mediante métodos numéricos, mediante integración por el método de los trapecios, 
pudiéndose expresar de la siguiente manera: 


La longitud 1 en la que existe el perfil, se obtiene conociendo un valor de h, y un 
incremento Ah , determinando el valor de h,. 
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Pavlovsky desarrolló un método para determinar el valor de / en canales, como 
resultado de la integración de la Ec. 97. Este método consiste en (Sandoval, 1993): 


1. Cuando , >0 


En donde: 
Na (e, -,) 
(1, =h,) 
K 
e = 2 
0 
£, 
ae, => 
£ 
Si h<h,: 
$=1,151l0g 2 
Si h> h,: 
+1 
$=1,151l0g2 
02= 
Ye: 
Pp? aC, ci B,, 
j SE 


2. Cuando 1 =(0 


En donde: 


8. Elsubíndice corresponde al valor medio de las características del flujo para las secciones 1 y 2. 
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En donde”: 
E (2, -e,) 
(», —h,) 
K 
e, = PR 
0 
K, 
e= 
K, 
), =arctane, 
Q, =arctane, 
Y: 
Há ac?-i-B 
Fr g:X 


El perfil de la superficie libre se determina conociendo por lo menos un punto el 
caudal y el tirante del flujo. Estos puntos se denominan “puntos de control”. La 
ubicación del punto de control depende del régimen del flujo. Así pues, para régimen 
sub-crítico las condiciones de aguas abajo son las que imponen la forma del reptil 
(remanso) y el punto de control. Para el régimen súper-crítico las condiciones de 
aguas arriba son las determinantes. 


9. Corresponde a un valor ficticio de , tomando el valor absoluto de la pendiente . 
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CAPÍTULO 3. PROBLEMAS 


1. Dado un canal trapecial con ancho de solera b=3 m, taludes 
laterales Mm, =1,5y M, = 2,0, pendiente longitudinal ¡ =0,0016 
y un coeficiente de rugosidad n= 0,013, determinar el gasto y la 
velocidad si el tirante normal h=2,5 m. 


Respuesta: O=71,26 m*/s; v=3,86 m/s 


2. Un canal trapezoidal de b=20m, 1,=0,0015, m=2 y 
n=0,025, transporta agua a razón de O = 40 mí / 5. Calcular el 
tirante normal y la velocidad del flujo. 


Respuesta: h=1,15 m; v=1,57 m/s 


3. Obtener la sección óptima y el costo por metro lineal (Cy de 
un canal trapezoidal de pendiente longitudinal ¡ =0,0010, 
coeficiente de rugosidad 1=0,015 y taludes laterales m, =1,5 
y Mm,=2,5, para transportar un caudal O=2 m 5/5. El canal ha 
ado revestido con un espesor de hormigón e=0,15 m, dejando 
un borde libre S =0,50 m entre la superficie de la lámina libre del 
agua y el nivel del terreno. Asumir un costo unitario del hormigón y 
excavación de 200 y 4 dólares por metro cúbico respectivamente. 


Respuesta: H =1,48 m ; B=1,01 m; C=241,75 dólares / m 


B 


4. Cuál es el caudal O que circula a través de una alcantarilla de 
hormigón simple (n=0,014) y de diámetro D=800 mm, si 
el tirante hidráulico en el tubo es d =0,30 m. Suponer que la 
alcantarilla ha sido instalada con una pendiente de 2 por mil. 


Respuesta: O =164,3 L/s 
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5. Obtener el diámetro de una alcantarilla de PVC (n=0,015) para 
transportar un caudal O =1 mis, instalada en una pendiente 
longitudinal ¿ = 0,001. Se desea que la alcantarilla funcione con 
un franco libre S =0,50 m sobre el nivel de agua. 


Respuesta: D=1,33 m 
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6. Enel canal trapezoidal que se observa en la Figura, calcular el ancho 
de la solera, tirante, perímetro y radio hidráulico, si la pendiente 
longitudinal del canal es 1 = 0,0001 y el caudal transportado es 
O=8 m? / s.. El canal se excava en tierra negra arenosa. Se sabe 
que la velocidad máxima permisible es v=0,5m/s ,elcoeficiente 
de rugosidad n= 0,02, el esfuerzo cortante T, =2,0 N/M yel 
ángulo de reposo O =357. 


Respuesta: b=11,52 m; h=1,19 m 
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